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Abstract. Applying Zeilberger's algorithm of creative telescoping to a family of cer- 
tain very-well-poised hypergeometric series involving linear forms in Catalan's constant 
with rational coefficients, we obtain a second-order difference equation for these forms 
and their coefficients. As a consequence we obtain a new way of fast calculation of Cata- 
lan's constant as well as a new continued-fraction expansion for it. Similar arguments 
can be put forward to indicate a second-order difference equation and a new continued 
fraction for £(4) = 7r 4 /90, and we announce corresponding results at the end of this 
paper. 



1. Introduction 

One of the most crucial and quite mysterious ingredients in Apery's proof [Ap], [Po] 
of the irrationality of £(2) and £(3) is the existence of the difference equations 



(n + l) 2 u n+1 - (lln 2 + lln + 3)u n - n 2 u n -i = 0, 



u 'o = !» u i = 3 > v 'o = °> v i = 5 > 



(1) 



and 



(n + l) s u n+1 - (2n + l)(17n 2 + 17n + 5)u n + n 3 u n = 0, 
u » = l, u '{ = 5, < = 0, < = 6, 
with the following properties of their solutions: 

lim 4 = C(2), hm % = C(3). 

"'n n 
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Unexpected inclusions u' n , D 2 v' n e Z and u", D^v" G Z, where _D n denotes the least 
common multiple of the numbers 1,2, ... , n (and Do = 1 for completeness), together 
with the prime number theorem {D]J n — > e as n — > oo) and Poincare's theorem, then 
yield the following asymptotic behaviour of the linear forms D 2 l u' n ((2) — D^v' n and 
Dn u nC(3) — D^v" with integral coefficients: 



<1, 



lim |^<C(2)-^<| 1/n = e 2 ( 
lim |^<C(3) - DWl 1/n = e 3 (V2 - l) 4 < 1, 

n— >oo 

and thus one obtains that both £(2) and £(3) cannot be rational. 

The two following decades after [Ap] were full of attempts to indicate the total list 
of the second-order recursions with integral solutions and to show their 'geometric' 
(i.e., Picard-Fuchs differential equations) origin [Za]. We do not pretend to be so 
heroic in this paper, and we apply quite elementary arguments to get new recurrence 
equations with 'almost-integral' solutions. 

In our general, joint with T. Rivoal, study [RZ] of arithmetic properties for values 
of Dirichlet's beta function 

1=0 v ' 

at positive integers s, we have discovered a construction of Q-linear forms in 1 and 
Catalan's constant 

i=o v ' 

similar to the one considered by Apery in his proof of the irrationality of C(2)- The 
analogy is far from proving the desired irrationality of G, but it allows to indicate the 
following second-order difference equation 

(2n + l) 2 (2n + 2) 2 p(n)u n+1 - q(n)u n - (2n - l) 2 (2n) 2 p(n + l)u n _! = 0, (2) 

where 

pin) = 20n 2 - 8n+ 1, 

(3) 

q(n) = 3520n 6 + 5632n 5 + 2064n 4 - 384n 3 - 156n 2 + 16n + 7, 



with the initial data 



7 13 

u Q = 1, u-i = -, v = 0, vi = — . (4) 
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Theorem 1. For each n = 0, 1, 2, . . . , the numbers u n and v n produced by the recur- 
sion (2), (4) are positive rationals satisfying the inclusions 

2 4n+3 D n u n E Z, 2 4n+3 D 3 2n _ lVn E Z, (5) 

and the following limit relation holds 

lim — = 67. 

The positivity and rationality of u n and v n follows immediately from (2)-(4). 
The characteristic polynomial A 2 — 11 A — 1 with zeros ( ( 1 ± y/E ) /2) J of the difference 
equation (2) coincides with the corresponding one for Apery's equation (1). Therefore 
application of Poincare's theorem (see also [Zul, Proposition 2]) yields the limit 
relations 

lim ul /n = lim vl /n = ( 1+ y^ l = exp(2.40605912 . . .), 

n^oo n^oo V 2 / 



lim \u n G - v n \ 1/n 



n^oo 



exp(-2.40605912...), 



2 

while the inclusions (5) and the prime number theorem imply that the linear forms 
2 4n + 3 £)| n _ 1 (it n 6? — v n ) with integral coefficients do not tend to as n — > oo. How- 
ever, the rational approximations v n /u n to Catalan's constant converge quite rapidly 
(for instance, |i>io/i*io — G| < 10~ 20 ) and one can use the recursion (2), (4) for fast 
evaluating G. Another consequence of Theorem 1 is a new continued-fraction expan- 
sion for Catalan's constant. Namely, considering v n /u n as convergents of a continued 
fraction for G and making the equivalent transform of the fraction [JT, Theorems 2.2 
and 2.6] we arrive at 

Theorem 2. The following expansion holds: 

13/2 1 l 4 -2 4 -p(0)p(2)| (2n-l) 4 (2n) 4 p(n-lMn + l)| 

[qi0) + \ ST) | qjn) 

where the polynomials p{n) and q(n) are given in (3). 

The multiple-integral representation for the linear forms u n G — v n similar to those 
obtained by F. Beukers in [Be, Formula (5)] for the linear forms u' n ((2) — v' n is given 
by 

Theorem 3. For each n = 0, 1, 2, ... , there holds the identity 

/• i /- i i"- i / 2 (i-i)» !/ "(i-i/)"- 1/2 
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2. Difference equation for Catalan's constant 

Consider the rational function 

7? m -nU<2t I n I ^ ~ I) ■ ■ ■ (t - n + 1) • (t + n + 1) ■ • • (t + 2n) 
R n (t) .- n\{2t + n + 1) {{t+ i ){t+ s y .. {t + n+ i ))3 

and the corresponding (very-well-poised) hypergeometric series 



t=o 



| r(3n + 2)r(n+i) 3 r(n + l) 



("I) 



r(2n+ f) 3 r(2n + l) 



X 6^5 



2'" 1 2 

Lemma 1. There holds the equality 



2 

Sn + 1, |n + §, ™+|, 71 +|' n + 1 

3 n + i 2n+ f , 2n + f , 2n + f , 2n + 1 



-1 



F n = U n (3(3) + U' n (3{2) + UMl) - V n 



where U n , D n U' n , D 2 n U'^ D\ n _ x V n G 2" 4 ™Z. 
Proof. We start with mentioning that 

p « 1)(t) := t(t -i)...( f - n + i) and 



^ 2) (*) == 



(t + n + 1) • • • (t + 2n) 



are integral- valued polynomials and, as it is known [Zu2, Lemma 7], 

2 2n -P n (-k- \) e Z for kez 



and, moreover, 



2 ^ D 3 . I djP n(t) 



eZ for keZ and j = l,2, 



t=-fc-l/2 



where P n (t) is any of the polynomials (10). 
The rational function 

Quit) 



(t+ !)(*+§). ..( t + n+ l) 



has also 'nice' arithmetic properties. Namely, 

(-l) fc (£)eZ if fc = 0,l,...,n, 
for other k G Z, 



a fc :=Q n (*)(t + fc + |)| 



t=-fc-l/2 
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that allow to write the following partial-fraction expansion: 

n 

1=0 1 + 1 + 2 

Hence, for j = 1, 2, . . . we obtain 

l-k 



f £(«.(«)(« + * + 5))l,._ t _i/2 = f £ X> (i - 



,=0 v t + i + 5 



t=-fc-l/2 



1=0 v y 

Therefore the inclusions (11)-(14) and the Leibniz rule for differentiating a product 
imply that the numbers 

A jk = A jk (n) : ~ (R n (t)(t + k + |) 3 ) \ t= _ k _ l/2 (15) 

= + n + !) • P n 1} (*) ' (*) • (Qn(t)(t + k + i)) 3 ) \ t= _ k _ 1/2 

satisfy the inclusions 

2 4n D 3 n ■ Aj k (n) E Z for fc = 0,l,...,ra and j = 0,1,2,.... (16) 

Mention now that the numbers (15) are coefficients in the partial-fraction expansion 
of the rational function (7), 

2 n . 

Substituting this expansion into the definition (8) of the quantity F n we obtain the 
desired representaion (9): 

2 j__ f_iy+k 



1= 0k=0 t=o^ t + ft+ 2) 



j=0 fc=0 
2 n 



j=0 fc=0 V Z=0 Z=0 7 v 2 

= t/ n /?(3) + C/;/?(2) + 0(l)-K 



where 



t/ n = 2 3 ^(-l) fc A 0fc (n), t/; = 2 2 ^(-l) fc A lfc (n), ^ = 2 ^(-l) fc A 2fc (n) 

fc=0 k=0 k=0 



(18) 



2 n k-1 , „w 



3=0 k=0 1=0 
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Finally, using the inclusions (16) and 

DL-iYl (oi , U3-i £Z for k = 0,l,...,n and j = 0,1,2, 

i =0 y u + L > 

we deduce that U n , D n U' n , D*U%, F)\ n _ x V n E 2~ 4n Z as required. 



Using Zeilberger's algorithm of creative telescoping [PWZ, Section 6] for the ratio- 
nal function (7), we obtain the certificate S n (t) := s n (t)R n (t), where 

s n (t) := — — — • (8n(2n - l) 2 (20n 2 + 32n + 13)t 4 

ny 1 2(2t + n + l)(t + 2n-l)(t + 2n) V v ;v ) 

+ 2(5440n 6 + 7104n 5 + 912n 4 - 1088n 3 + 76n 2 + 68n + 7)t 3 
+ (44800n 7 + 65600n 6 + 17568n 5 - 7056n 4 - 1088n 3 + 372n 2 + 146n - l)t 2 
+ (2n + l)(34880n 7 + 39328n 6 - 2176n 5 - 8416n 4 + 964n 3 + 154n 2 + 58n - 13)* 
+ 72(2n - 1)(2?2 + l) 2 (4720n 5 + 6192n 4 + 816n 3 - 864r2 2 + 69n + 13)) (20) 
satisfying the following property. 

Lemma 2. For each n = 1,2, ... , there holds the identity 

(2n + l) 2 (2n + 2) 2 p(n)R n+1 (t) - q(n)R n (t) - (2n - l) 2 (2n) 2 p(n + l)# n _i(t) 

= -S n (t + l)-S n (t), (21) 

where the polynomials p(n) and q(n) are given in (3). 

Proof. Divide both sides of (21) by R n (t) and verify the identity 

, 2/ v 2 / n / ,(2t + ?2 + 2)(t-?2)(t + 2?2 + l)(t + 2?2 + 2) . , 

(2n + l) 2 (2n + 2) 2 p(n) ■ (n + 1) 1 x1 ,J , 3v » " ^ 

(2t + n + l)(t + n + l){t + n + § ) 6 

— (2n — 1) (2n) p(n + 1) • 



n(2t + n + l)(t -n+ l)(t + 2n-l)(t + 2n) 

, ^ (2t + 72 + 3)(t+l) 3 (t + l)(t + 272 + 1) 

~ Sn[ + ) (2t + n + l)(t-n + l)(t + n + l)(t + n+lr ^ 
where s n (t) is given in (20). 

Lemma 3. The quantity (8) satisfies the difference equation (2) for n = 1,2, . . . . 

Proof. Multiplying both sides of the equality (21) by (—1)* and summing the result 
over t = 0, 1, 2, . . . we obtain 

(2n + l) 2 (2n + 2) 2 p(n)F n+1 - b(n)F n - (2n - l) 2 (2n) 2 p(n + l)F n _i = -S n (0). 

It remains to note that, for n ^ 1, both functions R n (t) and S n (t) = s n (t)R n (t) 
have zero at t = 0. Thus S n (0) = for n = 1, 2, . . . and we obtain the desired 
recurrence (2) for the quantity (8). 
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Lemma 4. The coefficients U n ,U' n ,U'^V n in the representation (9) satisfy the dif- 
ference equation (2) for n = 1, 2, . . . . 

Proof. We can write down the partial- fraction expansion (17) in the form 

4 +oo 



Rn{t) = E E 



A jk (n) 



j = l k= — oo v ^ ' 

where the formulae (15) remain true for all k G Z (not for = 0,1,..., n only). 
Now, multiply both sides of (21) by (—l) k (t + k + ^) 3 , take the jth derivative, where 
j = 0, 1,2, substitute t = —k — \ in the result, and sum over all integers k; this 
procedure implies that the numbers 

+00 +00 +00 

U n = 8 E (-l)*4>*(n), t£ = 4 E (-l) fc A lfc (n), 1^ = 2 £ (-l)*4ifc(n) 

k= — oo fc= — oo fc= — oo 

(cf. (18)) satisfy the difference equation (2). Finally, the sequence 

V n = U n p(S) + U' n (3{<2) + O(l) " F n 

also satisfies the recursion (2) by Lemma 3 and the above. 
Since 



2 3/4 3/4 7/4 7/4 

WW' {t) = " WW " + " WW- 



in accordance with (18), (19) we obtain 

U' = 8, U = U(f = V = 0, and ^ = 14, Vi = 13, ^ = U[' = 0, 

hence as a consequence of Lemma 4 we deduce that L/ n = U% = for n = 0, 1, 2, . . . . 
Lemma 5. There holds the equality 

Fn = U' n G — V n , 

where 2 An D n U' n G Z and 2 4n L>f n _ 1 V n G Z. 

The sequences u n := U^/8 and u n := V n /8 satisfy the difference equation (2) 
and initial conditions (4); the fact that F n ^ and \F n \ — > as n — > oo follows 
from Theorem 3 and asymptotics of the multiple integral (6) proved in [Be]. This 
completes the proof of Theorem 1. 
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3. Connection with 3 F 2 -hypergeometric series 

The corresponding very-well-poised hypergeometric series (8) at z = — 1 can be 
reduced to a simpler series with the help of Bailey's identity [Ba, Section 4.4, for- 
mula (2)]: 



1 + a — b — c, d, e 

3F2 , 

1 + a — b, 1 + a — c 



1 



r(l + a)r(l + a-d-e) 
T(l + a-d)T(l + a-e) 



x &F 5 



a, 1 + |a, 6, c, d, e 

|a, 1 + a — 6, 1 + a — c, 1 + a — d, 1 + a — e 



if Re(l + a — d — e) > 0. Namely, in the case a = 3n + 1, 6 = c = d = n+|, and 
e = n + 1, we obtain 



where the Euler-type integral representation for the 3F2-series can be derived as in 
[SI, Section 4.1] and [BR, proof of Lemma 2]. 
This completes the proof of Theorem 3. 

4. Conclusion remarks 

The conclusion (5) of Theorem 1 is far from being precise; in fact, using (2), (4) 
one gets experimentally (up to n = 1000, say) the stronger inclusions 1 

2 4n u n e z, 2 4n D 2 2n _ lVn e z. 

Unfortunately, they also give no chance to prove that Catalan's constant is irrational 
since linear forms 2 4n D\ n _ x {u n G — v n ) do not tend to as n — > oo. 

In the same vein, using another very-well-poised series of hypergeometric type 

and the arguments of Section 2, we deduce the difference equation 

(n + l) 5 u n+1 - r(n)u n - 3n 3 (3n - l)(3n + l)u n _i = 0, (22) 



1 A slightly weakened form of the inclusions is proved in [Zu4]. 
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where 

r(n) = 3(2n + l)(3n 2 + 3n + l)(15n 2 + 15n + 4) 

= 270n 5 + 675?i 4 + 702n 3 + 378n 2 + 105n + 12, (23) 

with the initial data 

uo = 1, u\ = 12, vq =0, = 13 
for its two independent solutions u n and v n , 

Theorem 4 [Zu3]. For each n = 0,1,2,..., the numbers u n and v n are positive 
rationals satisfying the inclusions 

6D n u n e Z, 6D 5 n v n e Z, (24) 

and there holds the limit relation 

i im ^ = L = ((4) = ri (25) 

n=l 

Remark. During the preparation of this article, we have known that the difference 
equation (22), in slightly different normalization, and the limit relation (25) without 
the inclusions (24) had been stated independently by H. Cohen and G. Rhin [Co] 
using Apery's 'acceleration de la convergence' approach, and by V. N. Sorokin [So] by 
means of certain explicit Hermite-Pade-type approximations. We underline that our 
approach differs from that of [Co] and [So]. 

Application of Poincare's theorem yields the asymptotic relations 

lim \u n \ 1/n = lim \v n \ 1/n = (3 + 2^3 ) 3 = exp(5. 59879212 ... ) 

n^oo n^oo 

and 

lim |u„C(4) -v n \ 1/n = |3 - 2 V / 3| 3 = exp(-2.30295525 . . . ), 

n^oo 

since the characteristic polynomial A 2 — 270A — 27 of the equation (22) has zeros 
135 ±78^= (3±2 v / 3) 3 . Thus, we can consider v n /u n as convergents of a continued 
fraction for £(4) and making the equivalent transform of the fraction we obtain 

Theorem 5. The following expansion holds: 

13 1 l 7 ■ 2 • 3 • 41 2 7 - 5 ■ 6 • 71 n 7 (3n - l)(3n)(3n + 1) I 

C 4 = HT^ + I TTT — + I — + • ' ' + I" ^ TT + • " ' 

|r(0) I r(l) !; r(2) \ r(n) 

where the polynomial r(n) is given in (23). 

The mystery of the £(4)-case consists in the fact that experimental calculations 
give us the better inclusions 

u n e z, D*v n e z 

(cf. (24)); unfortunately, the linear forms D 4 (u n ((4:) —v n ) do not tend to as n — > oo. 
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